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Abstract 

The resummation of 0(0^+^ In^' a;) terms in the evolution kernels of non-singlet combi- 
nations of unpolarized and polarized structure functions is investigated. The agreement 
with complete calculations up to order is demonstrated, and the leading small- a; contri- 
butions to the three-loop non-singlet splitting functions are derived. The additional 
contributions due to the resummed terms are studied numerically for the most important 
non-singlet structure functions. They are found to be about 1 % or smaller in the kine- 
matical regions accessible at present and in the foreseeable future. 
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1 Introduction 



The resummation of leading contributions in the evolution kernels of singlet structure function 
combinations at small x [1] may lead to large effects [2]. In this case, the small- a; evolution is 
dominated by the rightmost singularity in the iV-moment plane ~ [asl[N — l\f and higher order 

corrections to it. Such terms are absent in the non-singlet kernels both for the unpolarized 
and polarized structure functions [3, 4], as well as in the singlet kernels in the polarized case [5, 6]. 
Here the most singular contributions behave like N[as/N^y . An all-order resummation of these 
terms for non-singlet structure functions has been worked out in ref. [7]. Very sizeable correc- 
tions due to this resummation have been claimed for both unpolarized and polarized structure 
functions recently [8]. In this way the small- a; behaviour of the structure function evolution, e.g. 
of xF^'^{x, Q^), F^{x, Q^) — F^{x, Q^), and gi{x, Q^) — g'^{x, Q^), may be considerably affected. 

So far the resummation [7] was compared with the results of complete calculations only in 
the universal term of order a^jN [7, 8]. After setting up our notation and recalling the standard 
NLO formulation in section 2, we will show in section 3 that the resummation [7] agrees with 
the known evolution kernels P^s(a;,as) in next to leading order (NLO) for < as well in 
the small- a; limit. The contributions oc In^ a; to the so far uncalculated 3-loop non-singlet 
MS splitting functions are then derived from the results of ref. [7]. In section 4, we perform 
a numerical analysis for the most important non-singlet structure functions and compare the 
effect of the new terms beyond next to leading order with the NLO results. 



2 Evolution in fixed— order perturbative QCD 



The evolution equation for the non-singlet combinations f^si^jQ^) of part on densities reads 



(1) 



dlnQ^ 

Here (/) denotes the Mellin convolution, and P^s(a;, a^) is specified below. In the following, all our 
expressions refer to the MS factorization and renormalization scheme, and we drop the subscript 
'NS' wherever the non-singlet character of the quantity under consideration is obvious from the 
superscript 'it'. The splitting function combinations P^[x,as) are given by 



P^{x, as) = Pqq{x, as) ± Pqq{x, as) = ^ a[~^'^ Pf^ (x) 

1=0 

with as = as{Q'^)/{4:Tr). We now restrict ourselves to the spacelike case, 
expansion coefficients Pi~[x) obey the sum rule 

dxPfix) = V Z eAT , 



/ 

^0 



(2) 

> 0. The 



(3) 



which is due to fermion number conservation and Weierstrass' theorem, since acts as an 
independent parameter. At present the splitting functions are known up to two-loop order [3, 4]. 
They read 

'1 + x^' 



Pqq{x,as) = 2asCF 



+ 4af 
Pqq{x,as) = 4af 



C'fPf{x) + -CfCgPg{x) + CFNfTnPj,^{x) 



+ 0{a',) 



G'f 



-CfCg 



Pa{x) + 0{al) 



(4) 
(5) 
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where Cp = {N^ - 1)/{2Nc),Ca = N^, and Tr = 1/2. Nf denotes the number of active 
flavours. The functions Pi[x), I = F,G, Nf, A were derived in refs. [4]. For a; — )• the leading 
contributions to P^[x,as) are 



P+_o(a;,a.) = 2asCF + 2alCl\n^x + 0{a^^) 

P'^oi^, as) = 2asCF + 2al [-3C|, + 2CfCg\ In x + 0(a^) 



(6) 



The (scheme dependent) parton densities f^{x,Q^) are no observables themselves beyond 
the leading order. Instead of their evolution equation (1), one can directly consider the evolution 
of the structure functions F^[x,Q'^) obtained by the convolution 



F^H^,Q') = ct{x,Q')®ft{x,Q') 
Here cf[x^Q^) denote the coefficient functions 

CO 

ct{x,Q') = 8{l-x) + J2^[ct^M 



(7) 



(8) 



corresponding to F^[x,Q^). The evolution equation for Ff^ can be rewritten as an equation in 
as{Q'^) rather than in using 



da.. 



dlnQ^ 



This leads to 



dFt{x,as) 



where the NLO evolution kernels are given by 



Kf{x,as) (g) F^^{x,as) 



Kf^{x,as) = Pm,o{x)as + 



P±(^)-^Pns,o(^)-/3o<i(^) 
Po 



(9) 



(10) 



(11) 



Note the obvious fact that the In^ x terms of eqs. (6) enter the evolution equation (10) only in 
combination with the coefficient f3o in eq. (9). 

3 Resummation of leading sm.al\-x terms 

The transformation to Mellin-iV space 

M [k^^oM] (N) = jjxx^'-^K^^.ix.a,) = -lr±_^„(iV,a,) (12) 
of the most singular part of the evolution kernels in all orders in was given in ref. [7]^: 







> 






N2 


N dN^"" 


(e^^/^^-i/[2iV|](^)) 



(13) 



^Note that there are a few misprints in eq. (4.7) of ref. [7]. 
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Here z = N/^/2N^s, and Dp[z) denotes the function of the parabolic cylinder [9]. We expand 
eqs. (13) into a power series in aJ/iV^'*~^, and transform the result to a;-space using 



M 



InM- 

X 



{N) 



One obtains 



2asCF + 2alCl In' x + -a^C^ In^ x + 0(a^ In^ x) 

3 



(14) 



(15) 



2a, + 2a? 



In' X + -a^Cf 
o 



2iV2 



+ C( 



In^ x) 



The expressions (15) agree with the corresponding result found for P^_^q[x, as), eqs. (6), in the 
complete NLO calculations of the non-singlet anomalous dimensions [4] in the most singular 
terms since 

- = ^ + (16) 

holds in SU{N,). 

Besides the terms due to the anomalous dimensions Pf^[x), also the coefficient functions 
cfi{x) contribute in the evolution equation (10). The latter quantities have been calculated to 
0{al) for the structure functions F2{x,Q^), xF^^i^x^Q^) and g^i^x^Q^) [10, 11]. Expanding the 
coefficient functions for a; — )• 0, one finds that^ 



Ci^i{x) oc In ^-^ 



(17) 
(18) 



Therefore the terms of C(a^) and C(a^) in eqs. (15) can be identified with the parts of the 
non-singlet anomalous dimensions proportional to as(as In' a;)', assuming the validity of the 
resummation performed in ref. [7]. These contributions to P^[x) read 



-Cpln X 



2, X- 



10 



— C|. + 4(7|.(7g 



CfC^ ] In^ X 



(19) 



The calculation of the complete kernels K^[x,as) in eq. (10) in higher orders in a^ requires to 
take into account also higher orders in the /3-function. However, like in the NLO evoluOCon 
equation of section 2, the leading small- a; terms in oc as(as In' a;)' do not occur together with 
factors containing the coefficients /3i|i>i. 

It should be stressed that the agreement of the NLO terms between eqs. (15) obtained from 
the resummation [7] and eqs. (6) holds for g' < only. For the time-like case g' > one has [4]^ 



P+_o(a;,a.) = 2a,Ci.-2afC|,ln'a; + 0(a^) 

P'^oi^, as) = 2asCF + 2al [-5C|, + 2CfCg\ In' x + 0(a^) 



(20) 



^Note that apparent terms oc l/x™',m= 1,2 cancel in the corresponding expressions of ref. [11]. 

■^The 0[as) coefficient functions for e+e" annihilation or the Drell-Yan process behave also at most oc ln(l/a;), 



see ref. [10]. 
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The difference between eqs. (6) and (20) is due to the violation of the Gribov-Lipatov relation 
in the In^ x term of the NLO splitting functions. 

Considerations similar to those of ref. [7] may be valid for the spacelike twist-2 singlet evolu- 
tion equations in the polarized case, where the matrix of the splitting functions up to two loops 
has been obtained in refs. [5, 6]. 



4 Numerical results 

In moment space the evolution equation (10) for the non-singlet structure functions reduces to 
an ordinary differential equation. Taking into account the resummed kernels (13), the solution 
reads 



F±(iV,ao) 



7NS,o(JV)/2/3o 



(21) 



exp 



with 



If"" 1 I 
Jo 



+ 



ao 



Wo 



7i^(^)-1^7NS,o(iV) + 2/3oQ,i(iV) 



7f(iV) = -2 / dxx^-'Pt{x) , c: 



Jo 



dxx^-'4{x) 



(22) 



and ao = as{Ql). Here ^^{N) denotes the two-loop anomalous dimension with the term 
subtracted, since this contribution is accounted for already in the exponential factor. Moreover, 



r±(iv,a,) 



(23) 



The well-known solution in NLO for the evolution of F^[N,as) can be recovered from (21) by 
expanding the exponential to order a^. 

In the case of the non-singlet '+'-combinations the remaining integral in (21) can be done 
analytically, whereas it has to be performed numerically for the '-'-combinations involving the 
parabolic cylinder function. The transformation of the solution back into a;-space finally affords 
one standard numerical integral in the complex iV-plane [12]. We have also expanded the func- 
tions r^(iV, as) in the coupling constant a^. We find that in the practical cases considered below, 
one gets more than 90% of the resummation effect from the first two terms of the expansion. 

As it stands, eq. (21) violates the fermion number conservation for the '-' non-singlet combi- 
nations. Here the conjecture is that the coefficient functions cfi[x) do not contain terms oc In^' x. 
For this no proof exists yet, however, we have verified this behaviour up to 2-loop order in sec- 
tion 3 for the coefficient functions of xFs, , and g^^ . Under this assumption fermion number 
conservation has to be restored for T~^q[N, as). We approach this problem in two ways numeri- 
cally. In a first set of calculations we subtract a corresponding term oc ^(1 — a;) from the splitting 
functions P~, eq. (2), in each order in (the numerical results are labelled by 'A' later). In 
iV-space this prescription leads to 



T-{N,as)^T-{N,as)-T-{l,as) 



(24) 



Another possibility (denoted by 'B' in the following) is the restoration of fermion number con- 
servation by subleading 1/N terms modifying T~ according to 



T-{N,as)^T-{N,as)-{l-N) 



(25) 
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The difference of the results obtained by these two procedures gives an indication on the degree of 
dominance of the leading terms included in the present resummation vs. uncalculated subleading 
contributions. Our two prescriptions for restoring fermion number conservation are analogous 
to the procedure in the second reference in [2] with respect to energy-momentum conservation 
in the unpolarized singlet case. 

Before we come to the quantities studied numerically, we have to specify the input parton 
densities for F±(iV,ao) in (21). We choose = 4 GeV^ and A^js(iV/ = 4) = 230 MeV. In the 
present study we use the same input densities and value of Kqcd for the NLO and the resummed 
calculations. Specifically, in the unpolarized case we take the non-singlet combinations from the 
MRS(A) global fit [13]. For later use we note that xu^i^x^Q"^ behaves ~ x^'^^ at small x. In the 
polarized case, we employ Am„ and A(Z„ at = 10 GeV^ from ref. [14] ^ These densit les were 
obtained from the DFLM [15] unpolarized valence input distributions by multiplication with 
a;-dependent factors, yielding xAu^, ~ x°'^^ at small x. We have evolved these distribution back 
in NLO from = 10 GeV^ to our input scale Ql = 4GeV^ 

We are now ready to present the resummation effects on the most important non-singlet 
combinations. In the unpolarized case, we consider the evolution of the '-'-combination 

^ [xF^^'^ix, Ql) + xF^'^^ix, QD] = {x, Ql) ® [xu, + xd,]{x, Ql) (26) 
for an isoscalar target iV, and of the '+'-quantity 

F2''{x,Ql) - Fr{x,Ql) = 4^{x,Ql) (X) ^ [xu, - xd, + 2{xd-xu)] {x,Ql) . (27) 
In the polarized case we investigate 

gn^, Ql) - gri^, Ql) = ^o) ® I (A«„ - Ad,) {x, Ql) . (28) 

All these quantities are related to sum rules of phenomenological interest. xFs and gi are 
involved in the Gross-Llewellyn- Smith and Bjorken sum rules, respectively, which are used for 
the determination of a^- The F2 difference leads to the Gottfried sum rule, which provides 
information on the isospin asymmetry of the light quark sea in the proton. In all these cases, 
experimental data have to be extrapolated towards small a;, and hence the small- a; (^^-evolution 
is of interest in principle. A '+'-combination in the polarized case is the interference structure 
function g^^ , but this quantity is hardly measureable and the effects of the resummation are 
even smaller than in g^^ . We would like to note that the effects due to resummation in the 
'-'-combinations do not contribute to the first moment due to fermion number conservation. 

In figure 1, the NLO results for xFs and the resummed corrections AxF^'^ = xF^'^ — xF^^^ 
are displayed. The resummation effects turn out to be about 1% or smaller over the whole x 
and Q^ range considered. The inclusion of subleading terms to restore the fermion number 
conservation (curves 'B') reduces the effect by a factor of three or more at small a;, showing that 
even at x as small as a; = 10~^ the resummed terms do not dominate. With respect to the size 
of the resummation effect, the situation is the same for F2^ — F2"' as shown in figure 2. For this 
non-singlet combination there is no constraint due to fermion number conservation. 

The evolution of the polarized non-singlet structure function gl^ — g^ is depicted in figure 
3. Here the effect of the resummed Kernel (13) is bigger and reaches about 15% at a; = 10~^ for 
the (5^-values considered in the figure. This enhancement is due to the flatter small- a; behaviour 

^We use xo = 0.75 in eq. (12) of ref. [14]. 
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of the input densities {Auy ~ x~°'^^ vs. Uy ~ a;~°'^^)^. Note, however, that in the ^f^-range 
accessible for polarized electron anJ proton scattering at HERA [16] the effect is again about 1% 
or less. The measurement of the difference — g'^ will have even larger statistical errors than 
those of g^ alone, so that the effect due to resummation will not be resolvable in practice. With 
respect to the dominance of the resummed terms, the situation is the same as in the unpolarized 
case discussed for XF3. 

5 Conclusions 

We have investigated the effect of the resummation of terms of order a.['^^ In^' x given in ref. [7] 
on the small- a; behaviour of non-singlet structure functions for deep inelastic (polarized) lepton 
scattering both off unpolarized and polarized targets. The comparison with terms obtained in 
the same order by complete NLO calculations shows the equivalence of both approaches in this 
limit up to order a^. Since the coefficient functions up to two-loop order for the non-singlet 
combinations considered contain only terms less singular in In a;, the contributions to oc a^ln^a; 
in the three-loop splitting functions P^[x,as) can be predicted on the basis of ref. [7]. 

The numerical analysis shows that the all-order resummation of the terms 0[a.['^^ In^' x) leads 
only to corrections on the level of 1% for non-singlet structure functions accessible experimentally 
at present or in the foreseeable future. Moreover, for the '-' non-singlet combinations fermion 
number conservation has to be obeyed. This can lead to a further reduction of the effect by a 
factor of about three ore more in the small- a; range. The large sensitivity of the results on the 
prescription to implement the fermion number conservation constraint indicates the importance 
of so far uncalculated subleading terms 0{a['^^ In^ ~^ x) down to the lowst values of x considered 
here, x = 10~^. If it would be possible to measure the combination g^ — g'^ at high precision 
down to a;-values of this order, an enhancement relative to the NLO result ranging up to about 
5 — 15 % were obtained, where again the spread in the correction very roughly accounts for yet 
unknown subleading terms. 
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Figure 1: The small- a; (^^-evolution of the non-singlet structure function xF^ = ^{xF^^ -\-xF^^) 
for an isoscalar target N in NLO and the corrections to these results due to the resummed kernels 
derived from ref. [7]. 'A' and 'B' denote the two prescriptions for implementing the fermion 
number conservation discussed in the text. 
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Figure 2: The same as in Fig. 1, but for the combination — F^^. For this '+'-combination, 
there is no constraint on the kernels from fermion number conservation. 
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Figure 3: The small- a; (^^-evolution of the non-singlet polarized structure-function difference 
g^^ — g^"' in NLO and with the resummed kernels taken into account. Again 'A' and 'B' denote 
the two prescriptions for implementing the fermion number conservation discussed in the text. 
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